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ABSTRACT 

Context. Future observations of CMB anisotropies will be able to probe high multipole regions of the angular power spectrum, corresponding 
to a resolution of a few arcminutes. Dust emission from merging haloes is one of the foregrounds that will affect such very small scales. 
Aims. We estimate the contribution to CMB angular fluctuations from objects that are bright in the sub-millimeter band due to intense star 
formation bursts following merging episodes. 

Methods. We base our approach on the Lacey-Cole merger model and on the Kennicutt relation which connects the star formation rate in 
galaxies with their infrared luminosity. We set the free parameters of the model in order to not exceed the SCUBA source counts, the Madau 
plot of star formation rate in the universe and COBE/FIRAS data on the intensity of the sub-millimeter cosmic background radiation. 
Results. We show that the angular power spectrum arising from the distribution of such star-forming haloes will be one of the most significant 
foregrounds in the high frequency channels of future CMB experiments, such as PLANCK, ACT and SPT. The correlation term, due to the 
clustering of multiple haloes at redshift z ~ 2 - 6, is dominant in the broad range of angular scales 200 S / ^ 3000. Poisson fluctuations due 
to bright sub-millimeter sources are more important at higher /, but since they are generated from the bright sources, such contribution could 
be strongly reduced if bright sources are excised from the sky maps. The contribution of the correlation term to the angular power spectrum 
depends strongly on the redshift evolution of the escape fraction of UV photons and the resulting temperature of the dust. The measurement of 
this signal will therefore give important information about the sub-millimeter emission and the escape fraction of UV photons from galaxies, 
in the early stage of their evolution. 
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1. Introduction 

Three very important CMB experiments are due to operate in 
the ne xt two years: the PLANCK surveyor spacecraft (Taubei 
2004 ) and two ground based ex periments, ACT (Kosowskj; 
20061) and SPT dRuhl et allbooil) . They will cover the [150 - 
350] GHz frequency range with a resolution of a few arcmin- 
utes (I - 10^ - lO'*) and with unprecedented sensitivity. A lot of 
work has been invested in the last years in the study of possible 
foreground sources for these relatively high CMB frequencies. 
The galactic dust rad iation is obviously the most well 



ies, observable with the thermal and kinetic SZ ef- 



known foreground source (ISchlegel et al.lll998t iTegmark et al 
I2OOQ1) on the large and intermediate angular scales rel- 
evant for CMB observations, but it is very important 
to also look for foreground sources that have signifi- 
cant power at small angular scales (I ~ 10^ - 10"*). 
Several well studied foregrounds, like clusters of galax- 
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fects (ISunv aev & Zeldovichll972ll 1980t Komatsu & Kitavamal 



1999 



les 



Hernandez-Monteagu do et al 
199a 



2006h . SCUBA galax 
bright in frared galaxies 



(LS cott & White 19 99D and . „ 

dSanders & Mirabel 1996nLagache et al J2005h are predicted to 
contribute significantly to the Ci measured by the instruments 
mentioned above. It is relevant that these foregrounds are not 
only considered as difficulties impeding the study of primor- 
dial angular fluctuations, but that they might also provide us 
with extremely useful information on the statistical properties 
of objects like clusters of galaxies and compact sub-millimeter 
source s, and about the history of enrichment by heavy ele- 
ments (Basu et al. 2004). There are several projec ts, including 



lip: 

APEX (.GustenetaL.,2006,) . ASTE (lKohiiofr2005b and CCAT 



(ISebring et al.ll2006l) . which also plan to observe the sky in the 
sub-millimeter band. 

In this work, we will study the consequences of a sim- 
ple model in which haloes become bright in sub-millimeter 
wavelengths due to intense star formation resulting from merg- 
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ers and the c o rrespo nding supply of fresh gas. We use the 



Lacev & Cold (Il993h approach, based on the extended Press- 



Schechter formulation, to obtain the space density of the haloes 
in the universe. Knowing the space density of merging objects 
as a function of final mass and of the progenitor mass, we then 
introduce the lifetimes of the bright star formation phase as 
derived from numerical simulations of merging gas-rich disk 
galaxies. Using the Kennicut t relation betwee n star formation 
rate and infrared luminosity ( Kennicutt 1998 ) and taking into 
account that the bulk of infrared luminosity of the galaxies 
with intense star formation peaks in the sub-millimeter band 
at A. ~ 100 yum, we obtain the luminosity function of the merg- 
ing sub-millimeter objects. The emission spectrum of optically 
thin dust is a strongly increasing function of frequency (close 
to v2+/^, with yS = 1.5. Hildebrand 1983,) . Therefore distant ob- 
jects are observed at a rest-frame frequency closer to the peak 
of their spectrum. At frequencies below 1000 GHz, there is 
strong K-correction and the flux from galaxies at z > 1 remains 



appro ximately constant with increasing redshift (IB lain et al 



2002) . 



Recent observations of quasars at very high redshift (z ~ 6) 
and of extremely massive galaxies at z > 4 show that the 
chemical abundance of the gas (including dust) in the most 
dense parts of star-formin g galaxies is close to or exceeds 
the solar value ( FanI 2006 ). Therefore it is natural to assume 
that the dust radiation spectrum does not change dramatically 
with the redshift of the merging objects. This is the weakest 
assumption of the paper, which might lead to some changes 
in the dust emission spectrum at very high redshift when the 
abundance of heavy elements in gas and dust will drop below 
10"^ - 10"^ of the solar value. Nevertheless, observations of 
the Magellanic Clouds indicate that absorption and emission of 
dust remains similar to that in the Milky Way, when the abun- 
dance of chemical elemen ts is changing by one order of mag- 
nitude (Li & Drainell2002h . 

Our simple approach permits us to obtain several predic- 
tions that can be compared with the experimental results: 

1. the star formation rate density i n the universe as function 
of redshift (i.e. the Madau plot, Madau et al. I ll996lll998t 
lilopkins 2004; Hopkins & Beacom 2006 ); 

2. the upper limit on the intensity of the infrared backgr ound 
set by COBE/FIRAS observations jpixsen et al.lll998h ; 

3. the sourc e counts of sub-mil limeter sources performed by 
SCUBA dCoppin et al.ll2006h ; 

4. the redshift distributions of the brightest SCUBA sources 
C Chapman et al. 2005 ) . 

This model permits us to explain all but two of the above 
experiments. It demonstrates a deficiency of the brightest 
sources in comparison with the SCUBA data, but much more 
importantly it shows that the redshift distribution of the bright- 
est objects is also significantly different. We are then forced 
to make the simplest conclusion that SCUBA is observing ob- 
jects of a different nature. In this paper we are looking for the 
answer to the following question: where is the family of the sub- 
millimeter sources associated with galaxy mergings and why it 
is not detected yet in the direct source counts? Forthcoming 
SPT and ACT deep surveys will be able to find an additional 



population of sub-millimeter sources, since at the sensitivity 
levels of these instruments, up to ~ 50% of all detected ob- 
jects should be merging galaxies. This population should have 
a different redshift distribution than the SCUBA objects. Much 
deeper surveys with ALMA could also be dominated by sub- 
millimeter objects connected with merging galaxies, depend- 
ing on the properties of SCUBA sources, whose exact physical 
nature is still a matter of debate. Ongoing observations with 
SPITZER and the future results of Herschel will give addi- 
tional clues about the merging population in the infrared ob- 
servations. 

The merging model is able to produce a rather good fit to 
the SCUBA source counts if the dust temperature is close to 
25 K, but this value is not in agreement with the numbers pre- 
sented in the SCUBA papers. Additional energy coming from 
relatively weak AGN inside the merging objects could increase 
the flux from individual sources and mimic the source counts 
curve of SCUBA. 

The results of our computations show that the dust emis- 
sion of the merging objects should be one of the key foreground 
sources for primordial fluctuations of CMB, in the very broad 
range of angular scales I ^ 10^-10"*. Our model demonstrates 
the importance of the correlation term in providing the fluc- 
tuations due to haloes at high redshift (z ^ 2-6) and that 
its contribution exceeds the predictions for the Poisson fluctu- 
ations in the number of star f orming haloes. In this w ay we 
confirm earlier conclusions of Haiman & Knox (l2000b . about 
the importance of the correlation signal, which were based on 
a different model of sub-millimeter sources and on a different 
statistical characterization of the population (number of haloes 
instead of number of mergers). 

The correlation term from the relatively small mass (10** - 
10'' Mq) but numerous merging objects gives a contribution 
that samples the typical comoving correlation length of these 
objects (^ lO/i 'Mpc). This signal probes the most overdense 
regions in the universe at high redshift, that eventually will give 
rise to clusters and superclusters of galaxies. The Poisson term 
is smaller than the correlation term up to Z ~ 3000, but be- 
comes more important at higher I. ACT and SPT will be so 
sensitive that they will be able to perform successful source 
counts in their high frequency spectral bands which will permit 
the removal of the brightest sources (including mainly SCUBA 
sources) from the resulting maps at lower frequency, decreas- 
ing the contribution of the Poisson noise to the power spectrum 
of the background. This procedure will be unable to affect sig- 
nificantly the correlation term, which arises mainly from the 
numerous low-flux sources. 

We model here the sub-millimeter emission arising after vi- 
olent starburst episodes, with large populations of young stars 
appearing on very short timescales (of the order of 10^ - 10** 
years). The level of the predicted correlated signal depends 
strongly on the assumptions about: 

1. the escape fraction of UV radiation from merging star- 
forming galaxies and its evolution with redshift; 

2. the properties of the dust emitting in the sub-millimeter 
band (its effective temperature at high redshift and the 
emissivity index /?); 



M. Righi et al.: The clustering of merging star-forming haloes: CMB foreground 



3 



3. the redshift dependence of heavy elements and dust abun- 
dance; 

4. the amount of gas inside merging haloes. 

All these parameters are interconnected. Observations will give 
some hints about them and we will discuss below the informa- 
tion that could be obtained by precise observations of the an- 
gular power spectrum at different frequencies. 

Throughout the paper we use the recent WMAP3 
cosmological para meters for standard ACDM cosmology 
dSpergel et alj|2007l) . 



Major merger of two disks 



2. Halo distribution model 

To characterize the halo p opulation we will use th e standard 
Press-Schechter approa ch JPress & Schechtei 1974 ) in the ex- 
tended formulation of Lacey & Cold (Il993h . which also 
eludes a prescription for the halo mergers. 



m- 



2.1. Star formation from halo mergers 

We model the star formation within the framework of the halo 
model and derive the star formation rate as the rate of the bary- 
onic mass that is accreted per un i t time into new haloes. 

According to iLacev & Cold (1 19931) . the merger rate of a 
halo of mass Mi with another of mass Mi in a final halo of 
mass M - Ml + M2 is given by 



dN„ 



dMdt 



{Ml M,t) 



2 1 
n t 



d In 6ci- 



dlnt 



dcr 



dM 



6a- 



o-iMY 



1 



[1 -cr(M)2/o-(Mi)2]3/2 



X exp 



1 



1 



2 \cr(M)2 o-(Mi)2 



(1) 



where 6^ is the critical density for spherical collapse and 
cr(M)^ is the varia nce of the linear dens ity field. 
As pointed out by iBenson et al.l ( l2005h . the merger model is 
consistent only if the ratio between the masses of the merging 
haloes is not larger than ~ 10 - 100. Therefore in the follow- 
ing analysis we will include only major mergers, i.e. such that 

0.1 <: M1/M2 $ 10. 

In each merging episode a given amount of gas is converted 
into stars. We must take into account the cooling time which, 
for a given temperature T and metallicity Z, can be written as 



tcool(T,Z) - 



3 Pgas kT 

2^A(r,Z)«2' 



(2) 



3 c-1) 



where A(T,Z) is the cooling function (in units of erg cm-' s 
and yu is the mean molecular weight of the gas. The temperature 
of the gas is assumed to be equal to the virial temperature. The 
cool ing functions are tabulated for different values of T and 
Z bv lSutherland & Dopital(ll993h . For our estimate of the cool- 
ing time we choose a value Z = Zq. observations of quasars and 
galaxies at z ~ 6 bv lFanI (120061) show that the metal content of 
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Fig. 1. The star formation rate as a function of time during a 
merger of two gas-rich massive disk galaxies of equal mass 
(M ^ 10'^ /i^'Mq) used for the computations in this paper and 
based on the results of numerical simulations, cited in the text. 
Different colours identify diff'erent phases of star formation ac- 
tivity. 

such objects is close to the solar value, therefore this approx- 
imation is valid at least to this redshift. For even more distant 
objects the value of the metallicity is expected to be lower, but 
we will see that such objects will not contribute significantly to 
the effects that we are aiming to model. The cooling time sets a 
limit for the star formation, since only haloes with fcooi ^ fHub 
can cool within the cosmological time. This excludes the very 
massive haloes and the haloes with a mass corresponding to a 
virial temperature lower than Tvir - 10"^ K. 

Let us now consider two haloes of mass Mi and M2, merg- 
ing to yield a halo of mass M = Mi + M2. In each episode a 
given amount of the halo mass is converted into stellar mass. 
We parametrize this in the following way 

M2 



1 

Ml = — J] Ml 

* ' M/2 



(3) 



The amount of stellar mass formed from the halo Mi is there- 
fore proportional to its mass and to the mass of the merging 
halo M2, if this is sufficiently massive. In the same way, for the 
second halo 



2 Ml 

Mi - 77 M2 . 

* Qn, ' M/2 



(4) 



The total stellar mass produced in this merging episode is then 
the sum 



, 2 D.b Ml ■ M2 
M^^M\+mI^A — ti — -. 



(5) 



The parameter 77 is the star formation efficiency. We take this 
to be 5%, in order to match the observations of the cosmic star 
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formation history, which will be presented in Section l431 
Using this simple parametrization we are able to derive the stel- 
lar mass produced in a merger process between two haloes of 
given mass. 

2.2. Lifetime of tine starburst phase 

The approach described above allows us to compute the stellar 
mass produced in each merger episode. Since we are interested 
in deriving the rate at which gas is converted into stars and then 
connecting it to the luminosity of the sources, we introduce a 
characteristic timescale for star formation. 

Several numerical simulations of merging disk galaxies 
("Mihos & Hernauist 1994 Il996l: ISpringel & Hernquisd 12005 : 



Luminosity function 



[Robertson et al...2006.) show that the star formation process that 
occurs during a merger consists of multiple starburst episodes. 
Tidal interactions between merging galaxies can trigger star 
formation in the two disks before the final coalescence, al- 
though this usually coincides with the strongest burst. An accu- 
rate analytical description of such a process is complicated and 
thus numerical techniques must be used to find a solution. In 
the following, we will introduce a phenomenological prescrip- 
tion for the star formation in the merging haloes, based on the 
results of the numerical simulations mentioned above. These 
works studied the merging of two equal mass gas-rich galaxies 
and predicted an evolution of the star formation rate with time 
characterized by several phases, as displayed in Figure[T] 

1 . a long underlying phase (red), with a low and constant star 
formation rate, which is independent of the major merging 
process and is the sum of the quiescent star formation rates 
in each disk prior to the encounter; 

2. a second, more active phase (green), which corresponds 
to the first close passage of one galaxy around the other, 
when the tidal deformation experienced by each disk cre- 
ates shocks and compression of the gas. This triggers a 
moderate burst of activity of about 200 Mo/yr, with a dura- 
tion of the order of 3 x 10** years; 

3. a strong starburst phase (blue), which occurs when the two 
galaxies finally merge, with a star formation rate higher 
than 6OOM0/yr over a very short timescale, lower than 
10** years. 

We will take this model to be representative of the star forma- 
tion process during mergers of massive gas-rich objects. The 
bulk of the objects that contribute to the merging activity are at 
rather high redshift {z ^ 0.3-0.5), where the fraction of spirals 
and irregulars is higher than in our vicinity, where ellipticals 
are already a significant part of the total. 

As we have seen in Section IZTl each merger between two 
haloes M\ and M2 produces a certain amount of stellar mass 
Mi,(M\,M,z)- Such stellar mass will be created in the two ac- 
tive phases identified in Figure [T] in green and blue, accord- 
ing to the parameters summarized in the same plot. The first 
phase does not involve the stellar mass produced by the merger 
and hence will not be included in the model. Nevertheless, in 
Section lSTI we will estimate the contribution of this stage to the 
correlation term of the angular power spectrum. 
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Fig. 2. The shape of the far-infrared luminosity function of the 
merging haloes at different redshifts. 



We therefore have M* = M^^* -1- M]^' , with M^,-"' = 0.40 
and M\ - 0.60 M*. The percentages of stellar mass formed in 
each event are based on an average of the results presented in 
the literature and mentioned at the beginning of this section. 
Introducing the active time of each phase, it is easy to de- 
rive the corresponding star formation rate 



(3) 



f(2) 



where t. 



'acl 
(2) 



(0 



/ = 2,3. 



(6) 



2.5 X 



,8 (,-1 



yr and ffj 



:5x Wh-^yv. 



2.3. Infrared luminosity of the haloes 

The star formation rate given in (|6]l is strongly correlated with 
the inf rared luminosit y of the source, through the Kennicutt re- 
lation (lKennicuttlll998l) 



[Mo/yr] 



1.71x10-"'Lir[Lo], 



(7) 



where IR indicates the [8 - 1000] yum band. The bulk of emis- 
sion in this band peaks at ab out 100 jum and aroun d 80% of it is 
due to low temperature dust (i Lagache et al]|2005b . Therefore in 
our estimate we can directly use this relation and apply it to ob- 
tain the luminosity due to low temperature dust in star-forming 
galaxies. 

The number density of such bright objects for each population 
can be derived with the aforementioned model as 



dNl 



(0 

bright 



dV 



dVdt 



■{M,z)-t\ 



.(0 



(8) 
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Effects of the temperature evolution 
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Fig. 3. The effect of the correction in Equation ( flST l on the spec- 
trum of star forming haloes at z = 10: the solid line is the spec- 
trum with dust greybody only, the dashed line is the spectrum 
when including the coiTection due to the CMB. In the bottom 
insert the ratio between the two spectra at differrent redshifts 
is displayed to stress the importance of the correction with in- 
creasing redshift. The top insert shows the dust and CMB tem- 
perature evolution according to Equation ( fTSl i. 



where -j^{M\, M,z) is the number density of mergers per 
unit time, which can be obtained weighting Equation ([Hi over 
the mass function of the merging haloes 



dVdt 



iM,z) 



dM- 



dn dNmeig 

' dM, dMdt 



(9) 



remembering that M - Mi + M2 and that depends on both 
the masses of the merging haloes. 

The total number of bright objects at a given redshift will then 
be obtained by summing over the populations 



(2) 

bright 



dK 



(3) 

bright 



(M*,z). 



(10) 



dV dV dV 

Equation (fTOl i is a differential number density of bright objects 
in the universe. This is a crucial quantity, especially if we re- 
place M* with the infrared luminosity by taking advantage of 
the Kennicutt relation. In this way we obtain a statistical deter- 
mination of the sources responsible for dust emission, i.e. a lu- 
minosity function for these objects, which is shown in Figure|2] 
The sharp cutoff at Ljr ^ IO^-'Lq is due to the cooling time, 
which prevents the formation of the most massive (most lumi- 
nous) objects. 

2.4. Angular fluctuations from the star-forming objects 

The angular distribution of the emission of these star-forming 
sources will be characterised by their spatial distribution at 



different cosmic epochs. These sources will be located in the 
highest density regions in the universe, and their spatial clus- 
tering properties are well determined for a given cosmological 
scenario. Since most of the contribution comes from sources 
located at high redshift, their clustering properties turn out 
to be particularly relevant at angular scales of a few arcmin- 
utes. The shot noise associated with the Poissonian character 
of the source distribution will be relevant only at smaller angu- 
lar scales. 

In AppendixlAlwe provide a detailed computation of the statis- 
tical properties of the angular intensity fluctuations generated 
by the population of star-forming objects at high redshift. We 
find that, when probing sufficiently large scales, the change of 
brightness can be written, to linear order in the perturbations, 
as an integral along the line of sight 



Aly - const + 



r 

Jo 

r 

Jo 



dr 



C 47T C 



dr a(r) 6n{r) 



iv(l+z[r]) 3, . r\i-s:2\ 

— a (r) + 0(6 ), 

4-n 



(11) 



where r is the comoving distance, a(r) the scale factor and n(r) 
the number density of sources. The first term in the rhs is con- 
stant and introduces no anisotropy. The first integral contains 
the Doppler contribution, due to the peculiar velocity v of the 
sources, and the Integrated Sachs Wolfe (ISW) blueshift, due 
to time-varying potential ip, that the emitted photons experience 
as they travel towards the observer These two contributions are 
found to be negligible when compared to the second integral, 
which provides the contribution introduced by the varying spa- 
tial distribution of sources. This term can be separated in the 
contribution from the so-called Poisson term, accounting for 
the Poissonian nature of the objects, and from the correlation 
term, accounting for dependence of the source number density 
on the environment density. 

3. A model for dust emission 

A significant fraction of the stellar radiation emitted by galax- 
ies is reprocessed by dust and reradiated in the infrared and 
sub-millimeter band. The dust in star-forming galaxies has a 
very complicated distribution in space and a broad range of 
temperatures, depending on the proximity of dust grains to 
the sources of ionizing, optical and short-wavelength infrared 
radiation. The emissivity of the dust depends strongly on its 
composition and on the size distribution of grains. We do not 
have enough observational data to construct a reliable detailed 
model connecting the sub-millimeter spectra of the galaxies at 
redshift z > 2 with the star formation inside them. Several ob- 
servations (Blain 1999;,BIain et al. 2002; Chapman et al. 20051) 
show that in the majority of the brightest extragalactic infrared 
objects and in particular in SCUBA galaxies, the observed 
spectrum in the [100-3000] //m band is well fitted by a single- 
temperature modified blackbody v''Bv(7'du,st) with a tempera- 
ture Tdust - 30 - 50 K and an emissivity index in the range 
1 < y6 < 2; there is little change in the properties of dust along 
the Hubble sequence, thus we can expect this fit to be valid in- 
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Fig. 4. Testing our model with the SCUBA results. Different colours refer to the populations generated by the two bursts of 
star formation, as shown in Figure [T] (green and blue). Both plots assume Tdust = 30 K (with the CMB correction discussed in 
Section [TT), /3 - 1.5 and a mass-dependent escape fraction, which will be discussed in Section |4. 2. II Le^: the in tegral counts 
predicted by our model (black solid line) at 353 GHz (850 yum), compared with th e compilation of Coppin et al., (2006) fo r 
SCUBA sources (dots). The vertical lines represent s the sensitivity limits of SCUBA ( Coppin et al.l l 2006t Knudsen et al. 20061) . 
SPT and ALMA (expected in ~ 3 hours integration: iRuhl et al.ll2004l : iTarenghi & Wilsonll2005h . Right: the differential counts at 
353 GHz. We see that the short active phase is the most important and dominates the source counts at the high-flux limit. At the 
low fluxes the contribution from the two bursts is pratically equal. 



dependently of the Hubble type. As it is ob vious fr om the dis- 
cussio ns in the papers of .Reach et al.l (Il995i) and Schlegel et al. 
(Il998h . a fit with /? significantly smaller than 2 is attributable 
to multi-temperature dust or a difference in chemical composi- 
tion. 

For galaxies at high-redshift having similar rates of dust 
heating, the increase in the temperature of the cosmic mi- 
crowave background will affect the dust temperature. Therefore 
if the dust temperature in an object at z = is Tdust, fo r a similar 
object at redshift z it would increase as (lBlainlll999h 

Wz) = K.r-^((i..)-''-i)]""^r (12) 

where Tcmh - 2.725 K is the temperature of the background 
today. 

We consider the observations by PLANCK, SPT and ACT 
at fixed frequencies, namely 147, 217, 274 and 353 GHz; the 
peak of the SED for the majority of the brightest SCUBA and 
infrared galaxies is around 3000 GHz. This means that the fre- 
quencies higher than the peak will begin to dominate these 
observing channels at redshifts z ^ 19, 13, 10 and 8, respec- 
tively, so the flux observed in a given band will correspond to 
the flux emitted at a frequency higher than the peak. This has 
two consequencies. First, the contribution from high redshift 
objects will begin to decrease in the highest frequency bands of 
PLANCK, ACT and SPT. But in moderate frequency channels 
(147 and 217 GHz) these experiments continue to be sensitive 
to the first merging galaxies in the universe. Furthermore, for 
the high frequency channels the approximation of hot single- 



temperature dust emission will not be accurate, since in obser- 
vations of the local bright objects the region of wavelengths 
shorter than 100 fj.m shows that the decline in flux has a power- 
law behaviour, rather than the exponential cutoff of the single- 
temperature model. This demonstrates that the dust in reality 
has a broad distribution of temperatures which makes the cut- 
off after the peak much less steep. 

3.1. Dust emission in tiie CMB tiiermal bath at high 
redshift 

The fit by a modified blackbody spectrum considered above 
is valid under the further assumption that the influence of the 
CMB radiation field on dust grains can be neglected. Such an 
approximation is good at low redshift, when the background 
temperature is low, but as the redshift increases the intensity 
of this radiation becomes non-negligible. Therefore it is neces- 
sary to introduce a correction in the assumed spectral energy 
distribution of infrared sources 

<x / [B,(rdust) - Bv(rcmb)] , (13) 

where the difference in square brackets represents the correc- 
tion due to the presence of the CMB thermal bath. In the low 
frequency limit {hv <g: kT^asi) the luminosity is proportional 
to the difference between the grain temperature and the 
temperature of the CMB. As a result, for the same optical and 
ultraviolet energy absorbed by the grain, its radiation in the 
high frequency band will increase due to the existence of the 
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CMB thermal bath and correspondingly will decrease in the 
low frequency part of the spectrum. This is plotted in Figure [3] 
for an object at z = 10 with infrared luminosity Ljr - 10'"Lq. 



4. Three observational tests 

We identified the merger formalism presented in Section |2] as 
the main process leading to star formation in galaxies. Within 
the framework provided by the standard ACMD Cosmology, 
this approach constitutes a model for the IR emission in the 
Universe that can be compared with three different observa- 
tional constraints, as we show below. 

The model is particularly sensitive to the two spectral pa- 
rameters, the emissivity index /3 and the temperature Tdust of 
dust, which convert the bolometric luminosity into the spec- 
tral luminosity and hence fix the distribution of the infrared 
sources. Another important quantity, which will be discussed 
in this section, is the escape fraction of ionizing photons, which 
determines the fraction of ultraviolet radiation absorbed by dust 
and hence sets the relation between the star formation rate and 
the infrared luminosity. We will see that the introduction of 
such a parameter is necessary to satisfy the limits on the inten- 
sity of the infrared background set by the COBE/FIRAS data. 
Several observational tests can be performed to calibrate the 
free parameters of the model. 

4.1. SCUBA source counts 

Our merging halo model provides the source counts for the sub- 
millimeter so urces. Co mparison of our counts with SCUBA 
observations d Coppin e t al. 2006) at 850 fim (353 GHz) shows 
that, in the region above the SCUBA sensitivity limit at 0.2 
mJy, the curve is rather sensitive to the dust emission spectral 
parameters. In the left panel of Figure H] we plot the source 
counts resulting from our model. We will show in Section |6T| 
the variation of the curves for different values of the temper- 
ature and emissivity index: the closest result to the SCUBA 
counts is obtained for Tdust = 30 K and p = 1.5. These spectral 
parameters are consistent with those found from the observa- 
tions both in the local (Dunne & Bales 2001) and in the high 
redshift universe jChapman et al.l 120051) . so we will use them 
hereafter. 

The right panel of Figure |4] shows the differential source 
counts predicted by the model. As in the plot of integral counts, 
we show in different colours the contributions of the two ac- 
tive phases of star formation, presented in Figure[T] From these 
plots, it is evident that the short burst phase (blue) dominates 
the number counts, in particular at the highest fluxes, since 
it involves a comparable amount of stellar mass over a much 
shorter characteristic time. The rapid decrease of the differen- 
tial counts curve in the range [10 - 40] mJy reflects the rapid 
change in the merging rate occuring at z ^ 0.7 and possibly has 
the same nature as the rapid cosmological evolution of AGN. 

As an additional test of these fits, we plot in Figure |6] the 
distribution of the sources with redshift. In the left pa nel we 
compare with the observations of lChapman et al. (l2005h . in the 



same flux range [2.4, 17.4] mJy at 353 GHz. The two distribu- 
tions are significantly different, in particular our model predicts 
a peak around z ~ 1.5 while optical identifications of SCUBA 
sources give a peak at z ~ 2.5. This is an intrinsic limitation of 
the merging models, which predicts the formation of the most 
massive (and thus brightest) objects at redshifts lower than the 
redshift range of SCUBA sources. The short burst phase (blue 
line in the plot) plays the most important role, since it allows 
the detection of bright sources at much greater distances. The 
right panel shows several redshift distributions at the same fre- 
quency, for different flux ranges. 

Note, however, that the source number counts predicted by 
our model are slightly biased downwards when compared with 
SCUBA low flux data inferred from lensing measurements. We 
have found that, for the small sky area surveyed by SCUBA 
and the typical angular correlation scales for sub-millimeter 
sources, the source counts should show a typical dispersion 
much greater than the Poissonian estimate. In particular, we 
have found that typical fluctuations of the order of 20-25% of 
the real underlying number source count should be given when 
trying to estimate this quantity in areas comparable to that cov- 
ered by SCUBA. 

Observations and source counts in the submillimeter band 



have been also obtain ed by the MAMB O dBertoldi et al 



2007k and BOLOCAM (" Laurent et al]|2005b teams at 1.2 mm 
(250 GHz) and 1 . 1 mm (273 GHz), respectively. We compared 
their data with our predictions in Figure |5] we rescaled the 
BOLOCAM points to 250 GHz using the greybody spectrum 
assumed in this paper. The total sky area covered by these two 
instruments is comparable to that of SCUBA. The data of these 
three experiments are all consistent with each other. 

Most of the present and future experiments have a limited 
angular resolution (around 1 arcmin), which will not permit the 
resolution of two merging galaxies. The characteristic physical 
separation of two objects at their first passage is of the order of 
100 kpc (coiTesponding to 12 arcsec at z ~ 2) and of 10 kpc 
during the final coalescence (1.2 arcsec). But the impressive 
sensitivity and angular resolution of ALMA, which is expected 
to reach a valu e lower than 1 arcsec in the frequency bands 
discussed here jTarenghi & Wilson 2005), should enable both 
the identification of merging galaxies as double or asymmet- 
ric sources and the distinction between the two phases of star 
formation activity. Furthermore, ALMA will be able to demon- 
strate the difference between the typical merging starburst ac- 
tivity and the strong nuclear brightening due to the presence of 
a hidden AGN. 



4.2. Intensity of the cosmic infrared bacl<ground and 
comparison witli COBE/FIRAS results 

Another test that can be performed to calibrate this model is 
to compute the mean brightness of the cosmic infrared back- 
ground (CIB) due to the point sources and to compare it with 
the observations at the same frequency. Fixsen et al., (.1998.) . on 
the basis of the observations with COBE/FIRAS, 




Fi g. 6. Left: The redshift distribution of the sources at 353 GHz according to our model (soHd line) for the same flux range 
of IChapman et alJ (120051) and compared with their sample in the insert. Green and blue curves represent, as in Figure HI the 
contribution of the two star formation phases: also in this case the short phase is dominating and can be detected at much higher 
redshift. Right: the redshift distribution of the sources at the same frequency for different flux ranges: the [0.3 - 10] mJy interval 
(long-dashed line) could be relevant for SPT observations, the sources in the range [10^^ - 10"^] mJy (dotted line) contributes 
30% of the sub-millimeter background intensity at 353 GHz. 
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Fig. 5. The integral counts predicted by the merging model at 
250 G Hz compared with the results of MAMB O (.Bertoldi et all 
2007h and BOLOCAM jLaurent et al.ll2005b at flie same fre- 
quency. The points for BOLOCAM, originally presented at 1 . 1 
mm (273 GHz), have been rescaled to 250 GHz. 

presented the following fit for the intensity of the far-infrared 
background in the spectral range [150-2400] GHz 



where By is the Planck function and vq corresponds to 100 /vm. 
In our model we estimate the contribution of merging galaxies 
to the infrared background integrating the source counts over 
the flux 



rClB 



— I Sy dSv 

J dSy 



(15) 



We predict a large number of sources even at very small fluxes 
(5 353 ~ 10"^ niJy) and the slope of the differential counts in the 
range [10^^ - 10"^] mJy is so steep that these sources do not 
contribute much to the infrared background intensity. The bulk 
of it comes from the sources with flux in the range [10^^ - 10] 
mJy. However this is strongly related to the assumed value of 
the escape fraction of photons which strongly diminishes the 
contribution of less massive objects. 

4.2.1 . Escape fraction 

There is a natural process that decreases the contribution of 
star-forming galaxies to the intensity of the sub-millimeter 
background, namely the escape fraction of ionizing photons. 
Dust inside the haloes absorbes UV and optical photons and 
then re-emits them in the infrared band. However, assuming 
that all the ionizing radiation produced inside the galaxies is 
absorbed by the medium (i.e. by dust), so that the escape frac- 
tion is close to zero, is not completely correct, expecially for 
less massive objects, for which the escape fraction could be 
greater than zero. 

We model the escape fraction as a function of the halo mass 



rCIB 



(1.3: 



: 0.4) X 10"^(v/vo)° '"**''-'^Bv(18.5 ± 1.2 K), (14) 



/esc(M) = 0.8e 



(16) 
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Fig. 7. The spectrum of the cosmic infrared background pre- 
dicted by our model with a mass-dependent escape frac- 
tion (blue), Equation (fT6] l, compared with the fit to th e 
COBE/FIRAS observations (black) by iFixsen et all (Il998l) . 
Equation (fT4l i: dashed and dotted lines are the +lcr and ±2a" 
uncertainty levels, respectively. Red, green and orange lines 
corresponds to three alternative parametrizations for the escape 
fraction, see Section l673] for details. 



where Ms is the scaling mass, allowing an escape fraction of 
up to 80% for small-mass objects. To satisfy the limits imposed 
by the intensity of the infrared background we require a value 
Ms = lO'^/j-'Mo. 

The idea of a mass-depende n t esca pe fra ction was used 
in the works of Ricotti & Shull (l2000l) and iKitavama et al 



( 120041) . which dealt with the reionization of the universe. Here 
we are more concerned with the presence of dust which ab- 
sorbs ultraviolet photons and then influences the calibration of 
the Kennicutt relation, which assumes that almost all photons 
produced by young stars are absorbed by dust and re-emitted 
in the infrared band. In this sense, we need to compare the in- 
frared and ultraviolet luminosities of observed galaxies to de- 
rive an estimate of the escape fraction in these objec ts. 



For example, very recent observations bv lBuat et al.l (120070 us- 
ing SPITZER and GALEX allow us to estimate an escape frac- 
tion of the order of 1 - 2% for bright galaxies in the local 
universe. A very tight correlation bet ween the Lrm / ^Firv ra- 
and metallicity has been found by ICortese et al.l (12006) in 
a sample of nearby normal star-forming galaxies: the infrared 
to ultraviolet ratio increases with metallicity, thus leading to a 
higher escap e fraction for lo w metallicity objects. Earlier ob- 
servations by Heckman et al.l (fl998) seem to confirm this trend 
even for starburst galaxies, with escape fractions of the order 
of 30 - 50% for very metal poor objects. The metallicity is also 



Fig. 8. The evolution of the star formation rate according to the 
merger model (solid line), compared with the compilation of 
Hopkins & BeacomI (l2006h of numerous observations at differ- 
ent wavelengths (dots). 



known to correlate ve ry w ell with s t ellar m ass, as shown by 
Tremonti et al] ( |2004 and lErb et al.1 ( l2006h for the local and 
the high-redshift universe, respectively. In the merger model 
presented here stellar mass inside every halo is directly related 
to the mass of the halo itself. Then it is natural to expect a 
higher metallicity in more massive systems and therefore also 
a lower escape fraction. These simple arguments demonstrate 
that modelling this quantity as a function of mass is, at least 
qualitatively, well justified. 

As mentioned before, the Kennicutt relation implicitly as- 
sumes that almost all the ultraviolet photons produced by 
young stars are absorbed and reradiated by dust in the infrared, 
i.e. that /esc = 0. When considering non-zero values of /esc, we 
need to introduce a correction in the Kennicutt relation 



M„ [Mo/yr] 



1.71 X 10 



-10 



l-/es 



(17) 



TIR indicates the [1 - 1000] yum band 



In this way, less massive objects are unable to retain all their 
ionizing radiation, dust absorption will be less efficient and this 
results in a lower infrared luminosity for a given star formation 
rate. We will return to this point in Section l631 discussing sev- 
eral alternative parametrizations and their effect on the CMB 
fluctuations. 

As shown in Figure |7] with such a correction, our prediction 
for the intensity of the cosmic infrared background is within 
the lo" limits of the fit in Equation (fT4l l. On the other hand, 
with an escape fraction independent of m ass and redshift and 
close to the 2% value of Buat et al. IdaooA our model predicts 
an intensity for the sub-millimeter background which exceeds 
the COBE/FIRAS value but remains within the 2cr eiTorbars. 
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About 30% of the background intensity at 353 GHz comes 
from sources in the flux range ^353 - [0.01 - 0.1] mJy. The 
redshift distribution of such sources is displayed by the dotted 
line in the right panel of Figure |6] and it is characterized by a 
broad peak around z ^ 2.5 - 3. 

4.3. The Madau plot for cosmic star formation history 

The very simple parametrization which we use to model the 
star formation process within the framework of the halo model 
can be verified computing the cosmic star f ormation density as 



Star-forming haloes 



a function of redshift (i.e. the Madau plot , Madau et al.ll 19961 



19981: iHopkinsllIool iHopkins & Beacomll2006l) and compar- 



ing it with the observations. 

We average the star formation rate over all haloes Mi which 
are merging with other haloes to produce a final mass M 

r r dN^^,„ dn 
PM)^ dMi dM^—M.(Mi,M,,,), (18) 



where -j^- is defined in Equation ([T]) and M*(Mi , M2, z) is the 
number of stars produced in each merger episode. 

The evolution of cosmic star formation rate density is 
compared with recent observations, obtained at different 
wavelengths, in Figure [8] It is clear that our model slightly 
overestimates the cosmic star formation rate, especially at 
higher redshift (z ^ 5), where our points are outside the obser- 
vational errorbars. This could be due again to the limitations 
of the merging model which, though very accurate at the lower 
redshifts, fails to reproduce the merger trees at high redshift. 

Using a compromise model that is consistent with the measure- 
ments of the cosmic star formation rate that does not exceed 
the SCUBA number counts and that is close to the Icr limit 
of the sub-millimeter background provided by COBE/FIRAS, 
we can make predictions for the fluctuations in the CMB con- 
nected with the foreground sub-millimeter emission from star 
forming galaxies. 

5. Results 

In this section we present the angular power spectrum resulting 
from the distribution of dusty sources in the universe. As shown 
in AppendixlAl such a power spectrum is the sum of two terms, 
the Poisson Cj and the correlation Cf 



<A/,(«i)A/,(«2)) = 2 (cf + Cf)Pi(ni ■ H2), 



47r 



(19) 



where I is the multipole index and Pi is the Legendre polyno- 
mial of order /. 

The Poisson contribution to the angular power s pectrum can be 
obtained from the source number counts (e.g. IScott & White 
1999h 



Cj(y) 



(20) 



where all the sources brighter of 5 cut are supposed to be re- 
moved. We present the results for 5cut = 100 mJy at 353 GHz, 
unless otherwise specified. 
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Fig. 9. The Poisson (blue dashed line) and correlation (blue 
solid line) term of dust emission inside star-forming haloes, 
compared with the primary signal (black), at the four fre- 
quencies of the future CMB experiments (PLANCK, SPT and 
ACT): 353, 274, 217 and 147 GHz. Sources brighter ttian 100 
mJy at 353 GHz have been removed in the computation of 
Poisson term. 

The power spectra generated according to Equation ( |20] | 
and by ( IA.22I ) for the correlation term are shown in Figure |9] at 
the four different frequencies of the future CMB experiments. 
The correlation term peaks in a broad range / = [200 - 3000]. 
We notice that the Poisson signal is smaller than the correlation 
one for I ^ 3000 and becomes dominant only at the very large 
multipoles. Such multipoles, corresponding to angular scales of 
the order of a arcminute, are within the resolution limits of ACT 
and SPT. This plot shows that even at 147 GHz the fluctuations 
due to the dust emission in merging galaxies are one of the most 
important detectable extragalactic foregrounds. At higher fre- 
quencies it becomes dominant and at 353 GHz it should exceed 
the contribution from higher peaks (beginning from the fourth). 
As is well known, at 217 GHz the signal of the SZ clusters is 
negligibly small and this makes it much easier to observe the 
contribution from star-forming haloes. 

Since Poisson power is generated mainly by very bright 
sources (5^1 mJy), future experiments should be able to de- 
tect such sources in the maps and eventually to remove them. In 
this way the correlation term can be detected also at very small 
scales, since it is slightly affected by rare bright sources and 
obtains the bulk of its power from the more abundant and less 
bright sources. We estimate this in Figure [TO] at the 345 GHz 
SPT channel. This instrument is expected to have a sensitiv- 
ity of about 0.3 mJy, and this should allow one to excise the 
brightest sources and hence increase the weight of the correla- 
tion term versus the Poissonian one. In the same plot we also 
compare the predictions of our model with the result that would 




Fig. 10. The angular power spectrum of star-forming haloes at 
345 GHz: the red line is the predicted level of the fluctuations 
after the removal of the sources brighter than ~ 0.3 mJy (corre- 
sponding to SPT sensitivity, Ruhl e t al. 2004): compared with 
the model (blue) the Poisson term is reduced by more than one 
order of magnitude and the correlation term, which is much 
less dependent on bright sources, dominates up to higher mul- 
tipoles. The green line is the predicted power spectrum con- 
sidering only the sources detectable by SCUBA (with lensing 
technique) at 353 GHz (i.e. with flux greater than 0.2 mJy): the 
contribution of such sources to the correlation term is practi- 
cally insignificant, while they make all the Poisson power 



be obtained if only the sources detectable by SCUBA were con- 
sidered. The Poisson term would be at the same level, since it is 
generated mainly by the bright objects. The correlation would 
be strongly reduced, since SCUBA is not able to detect sources 
with flux lower than 0.2 mJy, missing the abundant dim sources 
which contribute significantly to the correlation power 

The analysis of the differential contribution to the signal 
from sources at different redshift reveals that the bulk of the 
correlation is generated by star forming haloes in the redshift 
range z ^ 2 - 4, with a slight dependence on the multipole I. 

Finally, following Ino ue & Kamava (|2004 . we estimate 
the amount of intergalactic dust expelled from galaxies by 
strong winds and AGN. This dust is heated by t he ultravi- 
olet b ackground and by the CMB. As shown by iBasu et aT 
(120041) . its effect on the CMB angular power spectrum will be 
6Ci ^ -2Tdust(v)C/ (high I), where C; is the primary signal and 
Tdust is the absorption optical depth. This turns out to be three 
orders of magnitude below the signature of the dusty merging 
haloes discussed in this paper 



Fig. 11. The e f fect o f foregrounds according to the model of 
Tegmark et al. I (l200flh . at 217 GHz: the long-dashed line refers 
to vibrating dust, free-free emission is plotted with the dotted- 
line and synchrotron radiation with the short-dashed line. The 
uncertainty on the noise level at 217 GHz for PLANCK is the 
red solid line, for a sky coverage /^j^y = 85%. Solid and dashed 
blue curves represent our results for star-forming haloes, the 
dotted curve is the correlation term from spiral galaxies (see 
SectionlSTTl. 

5.1. Contribution from normal spiral galaxies 

Another contribution comes from the infrared emission of nor- 
mal spiral galaxies, like the Milky Way. We decribe the distri- 
bution of such objects in space with the Press-Schechter mass 
function, weighted with the halo occupation number which 
gives the number of (bl ue) galaxies in a halo of given mass. 
Sheth & Diaferiol ( 2001 ) propose the following form for such a 
quantity 

A^gai(M) = J^j +0.5exp{-4[logio(M/10"-^5)]'}, (21) 

where M is the mass of the halo and with Mgai = 7x10'^/;"' Mq 
and a - 0.9. This fit is calibrated on the results of semi- 
analytical models and does not take into account the depen- 
dence on redshift. 

The correlation term resulting from this simple estimate is 
shown in Figure [TT] (dotted curve) and turns out to be more 
than two orders of magnitude smaller than the signal due to the 
star-forming haloes. 



5.2. Comparison with galactic foregrounds sources 

In Figure [TTI we present the well-known model of 
Tegmark etal J ( 2000l) for the main galactic foregrounds 
at 217 GHz (since at this frequency the contribution from 




Fig. 12. Sensitivity of the results for the integral source counts at 353 GHz to the choice of the spectral parameters: emissivity 
index {left, for fixed value of the dust temperature, 30 K) and dust temperature {right, for fixed value of the emissivity index, 
1.5). Grey dots are the SCUBA observations, as in Figure H] The curves in blue represent the choice which we use for all the 



results presented in t he paper. The plausible range of dust temperatures derived from the observations (iDunne & EalesI 12001 
Chapman et al]|2005h is 30-50 K. 
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Fig. 13. Sensitivity of the Poisson and correlation term at 217 GHz to the choice of the free parameters for the dust spectrum: 
emissivity index [3 {left, for fixed value of the dust temperature, 30 K) and dust temperature Tdust {right, for fixed value of the 
emissivity index, 1.5). The curves in blue represent the choice that we use for all the results presented in the paper. 



the thermal SZ is zero): the most intensive signal comes 
from the dust emission in our galaxy (long-dashed line), but 
it decreases rapidly at larger /, where our effect becomes 
dominant. The contribution from free-free emission (dotted 
line) and synchrotron radiation (short-dashed line) have 
similar dependence on I, but with much less power and should 
not affect the detection of the correlation term of dust in 



star-forming regions. 

We also include the expected uncertainty of the noise level 
for the PLANCK HFI detectors. At 217 GHz the expected 
resolution is 5 arcmin, with a sensitivity AT jT - 4.8 x 10"^. 
ACT and SPT telescopes are expected to give even better 
values of both resolution and sensitivity and hence a much 



M. Righi et al.: The clustering of merging star-forming haloes: CMB foreground 



13 



lower level of noise. This should allow the detection of the 
dusty sources in the power spectrum. 

5.3. Comparison with previous works 

An estimate of the expected angular p ower spec trum of the cos- 
mic infrared background was given bvl Haiman & Knox (2000). 
The model that they were using is different from the one pre- 
sented here: it assumes a Press-Schechter distributed popula- 
tion of sources with a mass and redshift dependent dust content 
and temperature related to the stellar field. Our model instead 
considers the distribu tion of merging objects, derived from the 
Lacey & Cold (Il993h approach, and thus directly probes the 
star forming haloes. Given these different approaches, we com- 
pare our result at 353 GHz with their "hot dust" model and find 
a rather good agreement: our correlation term at I - 1000 is a 
factor ~ 1 .4 (in Ar) lower than their estimate. 



Frequency dependence 



Song et al.l (l2003h performed detailed computations for 



the FIR background, modelling the dust emission from star- 
forming galaxies within the framework of the halo model, but 
not considering merging. They presented the angular power 
spectra for the 545 GHz channel of PLANCK assuming a 
sky coverage corresponding to the cleanest 10% and includ- 
ing the unsubtracted Galactic dust. Their results are essentially 

in agreement with Haiman & Knox ( 2000). 

More recently, lOonzalez-Nuevo et all (2005) presented 
their predictions for the angular power spectrum of extra- 
galactic point sourc es, based on the theoretical model of 
Granato et al.l (120041) . Comparing with their result at 353 GHz 
we find that at / = 1000 their correlation and Poisson terms 
are higher by a factor ~ 2.8 and ~ 1.6, respectively, than our 
predictions. 



Toffolatti et al.l (120051) presented an estimate of the power 



spectrum of ACBAR sources at 150 GHz. Our estimate for the 
Poisson term at the same frequency (and for the same flux limit 
■^cuf ~ 24 mJy) is roughly a factor 1.6 lower; we have a similar 
deficit when comparing with SCUBA observations. 

From these comparisons, it is clear that the contribution to 
the background fluctuations from the sub-millimeter sources 
in our model is lower than predicted by other authors, but in 
general their results are confirmed. 



6. Sensitivity to the parameters: can we learn 
something about the physical nature of the 
sources? 

The model presented in this paper is sensitive to several pa- 
rameters. As we pointed out in Section |4] it may be feasible to 
calibrate these parameters with the existing observations. The 
main purpose of this work is to make predictions for the level 
of foreground of dusty sources in the universe and to show that 
this signal will be important for upcoming experiments. 
However, this component should open an additional way to 
study the population of the most distant star-forming galaxies, 
which are already enriched by dust. The amplitude of the sig- 
nal could provide a very useful tool to understand the phys- 
ical nature of the objects. For example, the amount of dust 
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Fig. 14. The ratio of the correlation signal at different frequen- 
cies to the signal at 100 GHz, for different values of the dust 
temperature. 



could be derived if estimates are obtained about the actual 
value of the escape fraction, since these two quantities are 
related. Multifrequency observations will permit one to sep- 
arate the contribution of distant dusty haloes from other back- 
ground/foreground sources due to their characteristic spectrum, 
and could lead to a more precise estimate of the spectral prop- 
erties of dust emission, in particular the temperature. The ex- 
tremely high resolution of ALMA will make it possible even 
to individually detect very distant objects and to provide hints 
about their dust emission, in particular if information about 
their redshift is also retrieved (by observing, for instance, the 
emission in the fine-structure lines of the most abundant ions 
which were detected by COBE/FIRAS in the sub-milli meter 
spectr um of our galaxy and which were discussed in .Basu et alj 
(120041) ). Measuring the spectrum of the additional fluctuations 
in the CMB at different multipoles may also set constraints on 
the star formation scenario and confirm that the star formation 
at high redshift begins in strongly overdense regions which will 
eventually evolve into clusters and superclusters of galaxies at 
z ~ 0.5 - 0.7 and that the first star-forming haloes very rapidly 
create a high abundance of metals and dust. 
In this section we will explore the parameter space and present 
several results for the angular power spectrum as a function of 
different parameters. 



6.1. Spectral parameters 

The shape of the correlation term is fixed once the merger 
model is established, since it determines the distribution of star- 
forming haloes as a function of their luminosity. However, the 
merger model is connected only to the total infrared luminos- 
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Fig. 15. Left: the dependence of the correlation (solid) and Poisson (dashed) term on different parametrization for the escape 
fraction with redshift {right) at 353 GHz: low and constant (red), low and redshift dependent (green), high and redshift dependent 
(orange); blue lines correspond to the mass-dependent escape fraction adopted in the compromise model and show constant 
redshift evolution for given mass. 



ity through the Kennicutt relation. To obtain the value of the 
anisotropy in a specific frequency channel, one has to assume 
some kind of spectral distribution for the sources. In Section[3] 
we assumed a greybody distribution depending on the emissiv- 
ity index fi and on the dust temperature Tdust- Based on the com- 
parison with the SCUBA source counts (Figure fT2ll. we chose 
13-1.5 and Tjust - 30 K. These values are in agreement with 
the observational constraints, which give a wide range of val- 
ues both for the emissivity index (1.0 < < 2.0) and the dust 
temperature (30 K < Tdust ^ 50 K). As shown in Figure[T3] the 
amplitudes of the correlation and Poisson term are sensitive to 
the choice of the spectral parameters and display a variation 
over one order of magnitude in the range of possible values. It 
is clear, from these plots, that the two spectral parameters show 
a degeneracy. Multifrequency observations may allow one to 
reduce the extent of such degeneracy, as shown in the follow- 
ing. 

6.2. Frequency dependence 

As explained in Section |3] the contribution from high red- 
shift objects in the highest frequency bands of future experi- 
ments is expected to be low, due to the shift at wavelengths 
shorter than the 100 i-im peak in the emission spectrum. On the 
other hand, moderate frequency channels will be very sensi- 
tive even to extremely distant objects, since their emission will 
occur at wavelengths longer than the peak. This is shown in 
Figure [14] where the correlation term, normalized to the sig- 
nal at 100 GHz, is compared for different frequencies: while 
this ratio stays almost constant at 147, 217 and 274 GHz, it 
decreases at 353 GHz, towards the higher multipoles, corre- 
sponding to the smallest angular scales and hence to the most 



distant objects, whose emissivity at these frequencies is far be- 
yond the peak, on the Wien tail of the spectrum. 
The level of such a decrease at higher frequencies is somewhat 
temperature dependent. Higher dust temperatures shift the po- 
sition of the peak in the spectrum towards higher frequencies 
and shift the position of the sources beyond the peak towards 
higher redshift. This effect could be useful to measure the dust 
temperature of extremely distant objects, breaking the degen- 
eracy between and Tdust pointed out above, since the value of 
/3 has no influence on the position of the peak in the spectrum. 

6.3. Escape fraction 



As we mentioned in Section l4.2.1l we introduce an escape frac- 
tion for the ionizing photons in order to satisfy the limits of 
COBE/FIRAS on the intensity of the infrared background. In 
the compromise model, which provides all the main results of 
this work, the escape fraction is modelled as a function of the 
halo mass (EquationfTSIl. Here we present some other possibil- 
ities and show how these would affect our results on the fluctu- 
ations. 

GALEX and SPITZER recently obtained several estimates for 
the esca pe fraction in nearby luminous galaxies (L/g ^ 10" L©, 
see e.g. Buat et al.ll2007h . The published data do not give pre- 
cise information about the value of the escape fraction at high 
redshift (z ~ 3 - 5) from where most of the contribution to the 
correlation term comes. Similarly, up to now there is no clear 
idea about the dependence of the escape fraction on galaxy 
mass. We present here three further parametrizations: a con- 
stant and low escape fraction of 2%, an escape fraction grow- 
ing with redshift from 2% at z = to 40% at z = 4 and another 
one from 2% to 97%. 
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The coiTelation and Poisson terms at 353 GHz, coiTespond- 
ing to these parametrizations are plotted in the left panel of 
Figure [15] the blue line is the compromise model with the 
mass dependent escape fraction, the red line is the constant 
/esc = 0.02: the correlation term in this case in much stronger 
because the luminosity of the sources is weakly limited by 
the correction. The two Poisson terms are, on the other hand, 
very similar, since in the mass-dependent model the most mas- 
sive (and bright) objects that build up the Poisson power are 
also weakly affected by the correction. The green and orange 
lines are the two models with redshift dependent escape frac- 
tions: the coiTelation term is strongly reduced in the case of 
fesciz = 4) = 0.97, while the Poisson term, which is generated 
mainly by low redshift sources, is less affected. 

6.4. Gas fraction 

An additional parameter of the model is the gas fraction in- 
side haloes at different redshift. This quantity establishes the 
amount of gas that is available at each epoch for the formation 
of new stars. As explained in Section ITTl so far we have as- 
sumed that all the baryons inside the haloes are in gas form. 
This approximation is good at high redshift, when the number 
of stars is small, but it breaks down as more and more baryons 
are locked into stars. A simple way to describe this effect is to 
write the gas fraction at each redshift as 



/gas(z) = 1 - 



P*(Z) 



(22) 



where pb is the density of cold baryons inside haloes (the 
integral in the denominator is limited to those masses which 
are able to cool) and p*(z) is the stellar density, obtained 
integrating Equation ( fTSl ) in time. In this way we decrease the 
contribution of massive objects and reduce the level of the 
angular fluctuations about 30% (Figure [T6]l. 



7. Conclusions 

We have presented a simple model to determine the space 
density distribution of haloes which are bright in the sub- 
millimeter wavelengths due to intense s tar formation activity , 
resulting from mergers. We have used the Lacev & Cole ( 19931) 
approach to compute the merger rate and introduced a charac- 
teristic time of star formation for each merger episode, derived 
from numerical simulations of merging gas-rich disk galaxies. 
Using the Kennicutt relation we were able to derive the infrared 
luminosity of each object from its star formation rate. 

We assumed a dust emission model based on recent ob- 
servations of bright infrared galaxies by SCUBA (Blain et al. 
2002h . depending on two free parameters: the emissivity in- 
dex yS and the dust equilibrium temperature. We correct the 
spectral energy distribution taking into account the presence 
of the CMB thermal bath which becomes non-negligible for 
high-redshift objects. 

With these ingredients, we built estimates for several ob- 
servables and we calibrate the model in order to obtain a good 



agreement with available data: the Madau plot for cosmic star 
formation history and the source counts measured by SCUBA. 
Furthermore, we are within the limits imposed by FIRAS ob- 
servations of the Cosmic Infrared Background. 

In spite of these promising results, our model can be im- 
proved in several aspects. In particular, the merger rate ap- 
proach is known not to be accurate when compared to A^-body 
simulations. Our results depend crucially on several parameters 
(the escape fraction of ionizing photons, the dust temperature 
and emissivity index) which might be constrained very pre- 
cisely through observations of the angular fluctations, at differ- 
ent frequencies and angular scales, by future observations. We 
are confirming that the dust emission in merging star-forming 
galaxies is a powerful foreground for the high frequency chan- 
nels of PLANCK, ACT and SPT 
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Appendix A: The angular correlation of dust 
emission generated in star forming haloes 

In this section w e will follow the approach given by 
[Hernandez -Monteagudo et al. ( 20061) . 



A.1. Model of measured intensity fluctuations 

When considering sources that are much smaller than the PSF 
of the observing experiment, it is convenient to explicitly in- 
troduce the finite size of the beam and the frequency response 
of the detector(s). In what follows, our observing experiment 
will be characterized by a Point Spread Function (PSF) B(n) 
in denoting a direction on the sky) and a frequency response 
function 0(v', v), with v the central (observing) frequency. For 
simplicity, we shall assume that both functions have a Gaussian 
shape, so that 

\ 0' 1 

B(n) = exp -— , (A.l) 

b - 
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(^(v', v) = exp 



(y' - vf 



(A.2) 



The integral over the PSF and the spectral response will be re- 
garded as a normalization factor N: 



N= J dnB(n)dv' (p(v',v)^ 

- ( V27rO-b)^ y/TjTO-y = Qbeam Av. 



(A.3) 



The brightness intensity fluctuations introduced by the dust 
emission in haloes can be then written as 



AL 



^ /-fLSS 



drdnB(n)dv' (f>(v',v) 
+,[,.]) W(xj - r) a^(r), 



(A.4) 



where r denotes the physical distance from the observer and 
the index j sweeps over all star forming haloes in the universe, 
Xj being their positions; jv(\+z[r]) denotes the emissivity at the 
center of the halo placed at Xj, and W(Xj - r) is the dust profile 
and accounts for the dilution of the emissivity as the distance to 
the halo center increases. Note that due to expansion-induced 
redshift, the experiment will be sensitive to the emissivity at a 
frequency v(l + z{r\). Finally, a{r) = 1/(1-1- z{r\) is the scale 
factor at the epoch given by r, and the cr' factor accounts for the 
dilution of the intensity as the universe expands. We can replace 
the sum over j by a volume integral over the star-forming halo 
density nix): 

1 

77, 



A/, 



drdnB(n)dv' (f)(v',v) 
J dx n(x) jv(i+z[r]) W(x - r) a^(r). 



(A.5) 



Next we integrate in n and v, in order to obtain the contribution 
of all sources within the volume probed by the PSF and the 
frequency response of the instrument. The vector sweeping this 
volume will be y: we use the relation dv' - (v/r)dy, and the 
spectral response becomes a function of y, with the maximum 
aty = r. 

With this in mind, the argument of the function W is now x-y, 
resulting in 

A/v = dr J dxn(x)^^^ r^dndy(f>(y,r)jy(i^,i,^) 

X W(x - y) a\r) = 



AT Jo \r^) 4n 



jdxn(x)Vix-r). 



(A.6) 



In this case, the window function V{x - r) corresponds to the 
cosmological volume sampled by the finite size of the PSF and 
the frequency response. The integral on x corresponds to a con- 
volution that allows us to rewrite the last equation as 

r dk 

X 

J (2n) 



■ nkVk e 



(A.7) 



The Fourier modes of the halo density and the window V are 
given by and Vk, respectively. Using our model for the PSF 
and the frequency response, the latter can be written, for each 
r, as 



(A.8) 



The Fourier modes of the halo density will depend on the halo 
clustering properties, and will be characterised below. We can 
introduce the expression for Vu into Equation ( IA.7b . and this 
will remove the dependence on the normalization factor TV 



A/. 



Jo 



"■lss iv(l+:[r]) 3 

dr : a in 



r dk 

J (2^ 



dk 



4jt 



like 



-ikr 



(A.9) 



Further, if we are interested in scales larger than the beam size, 
the latter expression can be further simplified to 



AA 



Jo 



dr a (r) n{r). 

An 



(A. 10) 



which is simply an integral of the halo luminosity weighted by 
the halo density along the line of sight. From this expression, it 
is easy to see that the correlation properties of A/y will depend 
on the coiTelation properties of the star forming haloes. 
There are, however, two more aspects that need to be taken 
into account, namely the peculiar velocities of the star forming 
regions and the gravitational redshift imprinted on the emitted 
photons on the way to the observer. The former can be observed 
by simply imtroducing a Doppler term (containing the peculiar 
velocity of the sources) in the integrand of Equation ( lA.lOl i. the 
latter by explicitely writing the time-varying potential fluctua- 
tions in such integral. 
Thus, we rewrite Equation ( lA.lQb as 



My 



+ 2 



dr 



— a (r)n(r) 

An 



dr' \ \ + 



■^v(I+z[r']) 3 



An 



a\r')n{r') 



(A.ll) 



In this new expression, the term 0/c accounts for the time 
derivative of the potential, whereas the dot product v ■ n pro- 
vides the component of the source peculiar velocity along the 
line of sight. If we now argue that we are interested in scales 
where perturbations are still in the linear regime, then we may 
retain only the first order of the potential, velocity and density 
fluctuations that give rise to A/yi 



A/y - const 

"■'LSS 

+ I dr 

'LSS 



^«(.)^aV) + 2^/n.) 
c An c 



dr 5n{r) 
4;r 



'v(l+z[r]) 3 



where the average dust intensity Iy{r) is defined as 



7^'(r) 



r 



iir'Lv(i+z[r'])a {r')n{r') 



(A. 12) 



(A.13) 
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This equation shows that, to the first order, source peculiar ve- 
locities and potential fluctuations will couple to the average 
source number density «(r), and the correlation properties of 
the measured intensities will be ruled by the correlation func- 
tions of each of the three quantities: gravitational potentials, 
peculiar velocities and source number overdensities 6n(r). We 
shall find, however, that the latter will be dominant over the 
first two. 

A.2. Correlation of the measured intensity fluctuations 

We first rewrite the change in brightness intensity in terms of 
an integral of quantities expressed with respect to comoving 
coordinates 



Aly 



const + 



Jo 

/" 

Jo 



dr 



c An c 



dra(r) 6n(r) 



-v(l+zW) ^3 



4-7T 



a\r)+0(6-'). 



(A. 14) 



The comoving number density absorbs the a{r)^ factor, and a 
new a(r) factor appears from the line element. Note that pecu- 
liar velocities are proper velocities. Next, we express the star- 
forming halo number density as an integral over a luminosity 
distribution function 



dL 



dn 



'V(1+Z[AJ) 



dL 



(A.15) 



Finally, we decompose dn/dLy([+~i^]) as an integral over the 
haloes of different masses contributing to the same luminosity 
interval 



dn C dn ^ 

— ^ \ dM — G{M, L,(i+,[, 



]))■ 



(A.16) 



In this equation, the function G(M, Ly(\+~[x])) provides the frac- 
tion of haloes present in the mass interval [M, M + dM} (given 
by the Press-Schechter mass function dn/dM) that have re- 
cently experienced a major merger, and hence have given rise to 
significant star forming activity (see Section |2]l. We then write 
the spatial correlation function of haloes in the mass range 
[M,M + dMlni,(M,x), 

dn dn 
<«h(Mi,j:i)«h(M2,j:2)> = —(Mi,Xi)—(M2,X2) 
dM dM 

+ ^iMuXi)dl(Xi - X2)6DiMi - M2) 
dM 

dn dn 



environment. 

Having this present, one easily finds that the angular correlation 
function of the intensity fluctuations can be written as 

<A/,(ni) A/,(n2)> = ^ i^' + ■ m), (A. 18) 

with P/(ni ■ n2) the Legendre polynomia of order / and where 
Cj and Cf are the Poissonian and the correlation terms of the 
/-th angular power spectrum multipole. Cf is given by 



2 r r 

= - J dkk^ j 



•^r LSS y^_£2 

dra(r) ji(kr)— 

4-7T 



with 

In practice, this expression is equivalent to 



(A.19) 



(A.20) 



(A.21) 



where dN/dSy is the angular number density of sources per 
spectral flux unit. Due to its simplicity, we preferred to use 
Equation ( |A.21t in our computations. 

The correlation term can be expressed in terms of a A;-space in- 
tegral of the initial scalar metri c power spectrum PM ) times a 
squared transfer function, as in lSeliak & Zaldarriaga 



k^dkP^ik)\Ai(k)\\ 



(A.22) 



with the transfer function Ai(k) given by 



r 

A,ik) = 

Jo 



drjiikr) 



dS2(r) Vkic 



dr 



^s2ir)'^J-^n-' 



(A.23) 



In this equation, ji(x) is the spherical Bessel function of order 
I, 6k is the A:-mode of the dark matter density contrast, v^^ is the 
A;-Fourier mode of the peculiar velocity field of the emitting 
sources and i/r and ^ are the two scalar potentials perturbing the 
FRW metric (in the conformal Newtonian gauge). The func- 
tions S lir) and 52('") are defined by 

S\{r) = J dLy(u^[r])dMG(M,Ly^u^ir])) 
dn a(r)Ly(i+^[r]) 



dM 



An 



■b{M,z[r]) 



(A.24) 



and 



X b{Muz[x{[)b(M2,z[x2\)Uix\-X2). (A.17) ^2(0 = ^ dLy^i+,[,]-)dMG(M,Ly^i+.ir])) 

dn a(r)Ly^i+-[r]) 



The symbol 6d stands for Dirac delta and b{M, z) is the mass 
and redshift dependent bias factor tha t relates the halo an d the 
matter linear correlation function (IMo & WhitellT996l) . The 
first term on the right hand side accounts for the Poissonian 
fluctuations in the number counts, whereas the second term ac- 
counts for the dependence of the halo number density on the 



dM 



An 



(A.25) 



respectively. 



